Interpretation of the large-deformation high spin bands in selected A = 158 — 168 nuclei 
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The high-spin rotational bands in 168 Hf and the triaxial bands in Lu nuclei are analyzed using 
the configuration-constrained Cranked Nilsson-Strutinsky (CNS) model. Special attention is given 
to the up-sloping extruder orbitals. The relative alignment between the bands which appear to 
correspond to triaxial shape is also considered, including the yrast ultra-high spin band in 158 Er. 
This comparison suggests that the latter band is formed from rotation around the intermediate axis. 
In addition, the standard approximations of the CNS approach are investigated, indicating that the 
errors which are introduced by the neglect of off-shell matrix elements and the cut-off at 9 oscillator 
shells (J\fmax = 8) are essentially negligible compared to other uncertainties. On the other hand, 
the full inclusion of the hexadecapole degree of freedom is more significant; for example it leads to 
a decrease of the total energy of ~ 500 keV in the TSD region of 168 Hf. 

PACS numbers: 27.70.+q, 21.10.Re, 23.20.Lv 
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I. INTRODUCTION 

The high-spin structure of deformed nuclei shows a va- 
riety of interesting phenomena caused by the interplay 
between collective and single-particle excitations. The 
region of nuclei with Z ~ 72 and N ~ 94 is particularly 
fascinating. Potential energy surface (PES) calculations, 
predict that these nuclei constitute a new region of exotic 
shapes [IH3 coexisting with normal prolate deformation 
(e2 ~ 0.23). At high spins these nuclei may assume sta- 
ble triaxial superdeformed (TSD) shapes characterized 
by different moments of inertia for each of the principal 
axes. These TSD minima, with deformation parameters 
(e 2 ,7) ~ (0.4, ±20°), are caused by large single-particle 
shell gaps associated with proton numbers Z = 71 and 
72, and neutron numbers TV = 94 and 96 Hl|. Experi- 
mentally, such rotational bands have been reported in Lu 
(Z = 71) isotopes 

An extensive search for TSD bands in Hf (Z = 72) iso- 
topes has also been carried out, and a number of stro ngly 
deformed bands have been observed in 17 °- 175 Hf @4l3; 
where bands in 170 Hf Q and 174 Hf [ll[ have been tenta- 
tively assigned as triaxial. On the other hand, the pre- 
dicted TSD bands in 164 Hf and 166 Hf have not been dis- 
covered. Indeed, according to the analysis in Ref. fToj j . 
all observed strongly deformed bands in 17 °- 175 Hf are 
most likely near prolate falling into two groups corre- 
sponding enhanced deformation (ED) shapes (deforma- 
tions enhanced with respect to the normal deformed nu- 
clear shapes) and superdeformed (SD) shapes. The ED 
bands with £2 ~ 0.3 are built on the proton 213/2/19/2 
configuration while the SD bands involve the 7rz 13 /2 (pro- 
ton) and fj'15/2 (neutron) orbitals. On the other hand, a 
high-spin band has been observed in 168 Hf [14Ul6l] which 
appears to correspond to triaxial shape with a deforma- 
tion which is considerably larger than that of the TSD 
bands in 161 - 167 Lu. 

The high-spin bands which have attracted most in- 
terest recently are however the so-called ultrahigh-spin 



bands which bypass the band-terminating states in 
i57-i58 Er [7J and neighboring nuclei [H^. These 
bands were first assumed to have a triaxial deformation 
similar to that of the TSD bands in Lu nuclei but re- 
cent lifetime measurements [2l| show that they are more 
collective and they are suggested to correspond to either 
a larger triaxial deformation or possibly a similar defor- 
mation as the Lu TSD bands but with rotation around 
the intermediate axis (7 < 0). In a recent study [22l |. 
it was concluded that these bands must correspond to a 
larger triaxial deformation because the 7 < minimum 
appears to be a saddle point if the rotation axis is al- 
lowed to change direction. In any case, it has turned out 
to be difficult to find a consistent interpretation within 
the standard CNS approach (23hl25j . This is one reason 
why it appears important to investigate if, within the 
CNS approach, it is possible to get a consistent interpre- 
tation of the unique large deformation TSD bands which 
have been observed in 168 Hf. In this context, we will also 
demonstrate that the smaller deformation TSD bands in 
Lu isotopes appears to get a ready interpretation in the 
CNS formalism, see also (26j . 

Partly because of the large deformation of the TSD 
band in 168 Hf, some approximations of the CNS ap- 
proach become somewhat questionable. Therefore, we 
have made some modifications in the formalism making 
it possible to investigate the importance to include more 
oscillator shells in the basis and to account for all matrix 
elements coupling the different A^-shells of the harmonic 
oscillator basis. Most important however is that, for the 
first time to our knowledge, a complete minimization in 
the three hexadecapole degrees of freedom has been car- 
ried out at a large triaxial deformation. 

The motivation for the present work is to study high- 
spin rotational bands in 168 Hf and investigate their prop- 
erties in order to understand their nature. As a back- 
ground, we will consider the TSD bands in the Lu iso- 
topes. The Er bands have already been analyzed in 
Refs. [13, \2M but we will conclude with some additional 
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comments. We do the calculations within the frame- 
work of the configuration-constrained Cranked Nilsson- 
Strutinsky (CNS) model p3l - l25l | and another motivation 
is to test and develop this formalism. The model and 
standard approximations are explained in sect. [TTJ A brief 
description of the structure of the observed TSD bands 
in Lu isotopes using the CNS formalism is presented in 
sect. HIIl Standard approximations of the CNS formalism 
are tested in sect. IIVA1 while a complete minimization in 
the hexadecapole space is carried out in sect. IIVBI The 
reference energy which is often subtracted when present- 
ing nuclear high-spin bands is discussed in sect. IIV CI 
Then we study the experimental and theoretical high- 
spin bands in 168 Hf in sects. IV Al and IV Bl In sect. IV C| 
we compare these theoretical and the experimental bands 
and find out which theoretical bands correspond to band 
1, band 3, the ED band and the TSD1 and TSD2 bands 
of 168 Hf. Finally, we present some new points of view for 
the yrast ultrahigh-spin 158 Er band in sect. I VII 



II. THE STANDARD CNS FORMALISM 

In the configuration-de pen dent Cranked Nilsson- 
Strutinsky (CNS) model [23j-|25|, the nucleons are moving 
independently of each other in a deformed and rotating 
mean-field generated by the nucleons themselves. The 
rotation or the effect of the rotation is treated as an ex- 
ternal potential. The mean-field Hamiltonian used to 
describe a nucleon in the rotating nucleus is the cranked 
modified oscillator Hamiltonian 1231 



+ V 4 (e 4 ,7) - uj x 



(i) 



In this Hamiltonian, the cranking term uij x is introduced 
to make the deformed potential rotate uniformly around 
a principal axis with the angular velocity uj. The index 
t in the orbital angular momentum operator £ t . denotes 
that it is defined in stretched coordinates [13, HH . For 
168 Hf, standard values [23[ are used for the single-particle 
parameters n and /_/,, which determine the strength of the 
l t ■ s and £| terms, while A = 150 parameters [29| are 
used for the 161 ~ 167 Lu and 158 Er. This is motivated by 
the fact that the A = 150 parameters have been fitted 
for nuclei with N sa 90, while standard parameters should 
be more appropriate for the well-deformed nuclei in the 
middle of the rare-earth region. 

In Eq. (p}, hiio{^2il) i s an anisotropic harmonic- 
oscillator Hamiltonian: 



P ,1 / 2 2 2 2 2 2, 

— + -771(^0; +oj y y +oj z z) (2) 



hHO (£2,7) 

The relation between the oscillator frequencies and £2 , 7 



u x = w (e 2 , 7) (1 ~ ^£2 cos( 7 + y )) 

ui y =w (e 2 ,7) (1- cos( 7 - y)) 
2 

w z = w (e 2 ,7) (1 - g£2 COS7) 



(3) 



The deformation dependence of O2o(£2,7) is determined 
from volume conservation of the equipotential surfaces. 

The total energy is obtained using the shell correction 
method. Thus the shell e nergy , E s h, is calculated using 
the Strutinsky procedure [30, [31 1 and the total energy is 
defined as the sum of the shell energy and the rotating 
liquid drop energy [H,[3l|, E rU , 



Etot(I) = E sh (I) + E rld (I). 



(4) 



This renormalization ensures that the total nuclear en- 
ergy is correct on the average. The Lublin Strasbourg 
drop model (32j is used for the static liquid drop energy 
with the rigid-body moment of inertia calculated with 
a radius parameter rn == 1.16 fm and a diffuseness pa- 
rameter a — 0.6 fm [251 ] - Finally, minimizing the total 
energy for a given angular momentum with respect to 
deformation gives the equilibrium shape and correspond- 
ing energy. Plots of the minimized total energy versus 
spin I are frequently used in the description of high-spin 
properties of rotating nuclei. To present considerably 
more detailed information about individual and relative 
properties of the rotational bands, the excitation energy 
is plotted relative to a reference energy. Note that the 
same reference energy is utilized for all theoretical and 
experimental energies in a nucleus. 

Eq. ([1]) represents the rotating modified oscillator 
Hamiltonian in terms of the quadrupole, £2, non-axial, 7, 
and the hexadecapole, £4, deformation parameters. The 
dependence of the Hamiltonian on the hexadecapole de- 
formation is written as: 



with [23l . [3 



y 4 = 2^ p 2 [£40*4 {6t,<Pt) 

+ e 42 (r 4 2 (# t ,^ t ) + y 4 - 2 (6 t ,tpt)) 
~*~ £44 (r 4 4 (0t,¥>t) + *7 4 (0*>^))] - 



£40 = £4^(5 cos 2 7 + 1) 
6 



£42 



(5) 



-£4— v 30 sin 27 



(6) 



1 , 

£44 = £4 — V 70 sin 7, 



where 9t and tpt are the polar and azimuthal angles in 
stretched coordinates and p is the radius in stretched co- 
ordinates. The 7 dependence in Eq. © is introduced in 
such a way that the axial symmetry is preserved when 
7 = -120°, -60°, 0° or 60°. All ellipsoidal shapes can 
be described within a 60° degree sector, but the rota- 
tion occurs around the shortest, the intermediate and the 
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longest principal axis for 7 = [0°, 60°], 7 = [0°, -60°] 
and 7 = [-60°, -120°], respectively. 

Because the en parameters depend on one parame- 
ter £4, there is only one hexadecapole degree of freedom. 
In a standard calculation, the total energy is minimized 
varying three parameters: two quadrupole parameters, 
e 2 and 7, and one hexadecapole parameter, £4 [23|. The 
choice of the deformation space to be used in a calcula- 
tion is important. Recently, some studies concentrating 
on the role of different multipoles on the fission barrier 
heights have considered more general hexadecapole de- 
formations [33-l36|. 

The rotating basis \n x n 2 n 3 Ti) can be utilized to diag- 
onalize the Hamiltonian matrix and to find eigcnfunc- 
tions of Eq. (fTJ) [23| ■ Since the couplings of j x are fully 
accounted for in the rotating basis, the only terms in 
Eq. (TTJ) which couple between basis states of different 
Nrot — n x + fT-2 + n.3 are the hexadecapole deformation 
potential V4, and the it ■ s and £ 2 terms. The off-shell ma- 
trix elements of the latter terms are small for reasonable 
rotational frequencies. The importance of the off-shell 
matrix elements of the V4 term depend on the deforma- 
tion region where hexadecapole deformations generally 
become more important with increasing quadrupole de- 
formation. For small £4 values it thus seems reasonable 
to neglect all those matrix elements which are off-shell 
in the rotating basis and keep N ro t as a preserved quan- 
tum number. The important advantage of the rotating 
basis is that N ro t (generally referred to as N below) can 
be treated as an exact quantum number making it pos- 
sible to fix configurations in great detail. It seems that 
this is the most important feature explaining the success 
of the CNS approach; especially the possibility to follow 
e.g. terminating bands in spin regions where they are not 
yrast. 

The diagonalization of the Hamiltonian, Eq. (flj , gives 
the eigenvalues , which are referred to as the single- 
particle energies in the rotating frame or the Routhians. 
Subsequently, it is straightforward to calculate different 
expectation values like (j x ) and (j 2 ). The diagonaliza- 
tion of the Hamiltonian is performed with a cut-off in 
the single-particle basis which may lead to errors in the 
results. The original CNS codes were written with only 
9 oscillator shells (J\f ma ,x = 8) in the basis and this is the 
maximum number of shells which has been used in all 
subsequent CNS calculations, e.g. [I?], H3, [53] ■ It seems 
important to test these approximations, i.e. the neglect 
the off-shell hexadecapole matrix elements and the cut- 
off in the rotating single-particle basis. 

In the present calculations, pairing correlations are ne- 
glected, although, it is quite evident that the pairing field 
is essential for the description of atomic nuclei 38] . This 
is seen for example from the observed energy gaps and 
the suppression of the moments of inertia in rotating nu- 
clei. However, it appears that the most of the properties 
of nuclei at high spins are rather insensitive to the pairing 
field. For example, rotational bands have been studied 
by the Cranked Nilsson-Strutinsky approach [23T - I25] . the 



Cranked relativistic mean field theory (39l - l4l| not includ- 
ing pair correlations and the Cranked relativistic Hartree- 
Bogoliubov formalism [42T444T ] including pair correlations. 
These studies show that in high-spin regime, calculations 
without pairing describe the data accurately. In view of 
this, it is often advantageous to carry out calculations in 
an unpaired formalism because of the more transparent 
description and, for the present CNS calculations, the 
unique possibilities to fix configurations, making it pos- 
sible to follow for example the drastic shape changes in 
terminating bands [24], [45[ . 

In order to evaluate the importance of the pairing 
energy in the odd-odd 76 Rb nucleus, rotational bands 
have been studied by the Cranked Nilsson-Strutinsky- 
Bogoliubov (CNSB) formalism presented in Ref. [46| with 
particle number projection and with energy minimiza- 
tion not only in the shape degrees of freedom, £2, 7 and 
£4 but also in the pairing degrees of freedom, A and A 
and have been compared with the predictions of the CNS 
model [47j . In these calculations, the contributions from 
pairing are found to be small at low spin values and they 
decrease with increasing spin. The pairing energies do 
not change the general structure which means that, for 
example, the potential energy surfaces with pairing in- 
cluded are found to be very similar to those in the CNS 
formalism. 

The outcome from CNS and CNSB calculations have 
also been compared in 161 Lu [H,[4§|. It turns out that for 
/ > 30, the inclusion of pairing will correspond to a small 
renormalization of the moment of inertia but it does not 
affect the general structure of the yrast line, band cross- 
ings etc. Especially, the terminating states for / ~ 50 
are essentially unaffected by pairing correlations. With 
this in mind, we will analyze the high-spin states of Lu 
isotopes and 168 Hf in the unpaired CNS formalism where 
our main interest are those configurations which cannot 
be isolated in present formalisms with pairing included. 

For A = 158 — 168 nuclei, it is convenient to la- 
bel the configurations by the dominant amplitudes of 
the occupied orbitals and holes relative to the 146 Gd 
(Z = 64, N = 82) closed core; that is, 

K(hn/2) Pl (h 9/2 f 7/2 yz(i 13/2 y 3 

V (M = 4 )-»* (/Mi/ 2 )-" 2 ( i i 3 /2)" 3 (in/2, 59/ 2 )" 4 O15/2)" 5 , 

where the number of the Af = 4 protons and h s / 2 , fr 12 
neutrons is determined from the total number of pro- 
tons and neutrons in a nucleus. We will often use the 
shorthand notation (where the numbers in parentheses 
are omitted when they are equal to zero), 

[P1CP2P3), (nin 2 )n 3 (n 4 n 5 )}. 

Note however that this is only for the purpose of labeling 
the configurations; in the numerical calculations no core 
is introduced and all or most of the couplings between 
j-shells are accounted for according to the different ap- 
proximation schemes. 
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III. TSD BANDS IN LU ISOTOPES 

The TSD bands in Lu nuclei are characterized by an 
odd H3/2 proton which plays an important role for the 
wobbling excitation Q. Apart from this, the occupied 
orbitals in these bands have not been given much atten- 
dance. An exception is Ref . || , where the single-particle 
orbitals and the corresponding shell gaps at TSD defor- 
mation were discussed. Here we will try to demonstrate 
the filling of the orbitals in the lowest TSD bands, in- 
dicating the contribution of the, specific orbitals which 
become occupied when the number of neutrons increases. 
This is analogous to previous classifications of the su- 
perdeformed bands in the A — 150 region [37], H3, [H| . A 
preliminary report of the present classification was given 
at the NS2008 conference [Hj]. 

As seen in Fig. 5(a) in Ref. [2(| (and in Fig. [TBI be- 
low drawn at a somewhat larger deformation), the pro- 
ton configuration with two hg/2 and one H3/2 proton is 
favoured for TSD deformations (£2 ~ 0.37, 7 ~ 20°) for 
frequencies up to fiw ~ 0.6 MeV. Indeed, according to 
our calculations, this is the proton configuration, 8(21), 
for the lowest calculated TSD bands in the 161_167 Lu 
isotopes. In order to understand the neutron configura- 
tions, Fig.[T]is instructive. Starting from the left, it shows 
the single-neutron orbitals for prolate shape in the range 




0.15 0.25 0.25 0.35 0.45 
£ 2 Y=0° 10° 20° e 2 



FIG. 1. (Color online) The single-neutron orbitals drawn 
along a path in the (e2,7)-plane in order to clarify the origin 
of valence orbitals in the TSD minimum of A = 160 — 170 nu- 
clei. The orbitals are labelled by the approximate asymptotic 
quantum numbers, fi[A/Vi z A], also for 7/0, even though 
fi, which is the projection of j, is not preserved in this case. 
The upsloping orbitals emerging from the subshells below the 
N = 82, which are important when building strongly col- 
lective bands, are highlighted. Note how the 1/2 [400] and 
3/2 [402] orbitals repel each other with increasing axial asym- 
metry, thus inducing triaxial shape in configurations with 
holes in the 3/2 [402] orbital. The N = 92 shell gap of the 
TSD band in 163 Lu is marked out and it is then shown how 
the TSD bands in 161 ' 185 > 16 'Lu are formed from holes or par- 
ticles in the valence orbitals. 



£2 = 0.09 — 0.25, then for £2 = 0.25 as a function of axial 
asymmetry 7 and finally for constant 7 = 20°, again as 
a function of £2. The neutron configurations of the TSD 
bands in the Lu isotopes with N = 90 — 96 are then illus- 
trated at £2 ~ 0.40 (and 7 = 20°). The gap indicated for 
N — 92 is responsible for the 163 Lu configuration which 
has two holes in M — 4 and two holes in hxi/2 Af = 5 or- 
bitals combined with six particles in AT — 6 orbitals, i.e. 
the configuration (22)6. As discussed e.g. in Ref. [l7j|. the 
holes in the upsloping Af = 4 and ft.11/2 orbitals are very 
important for the formation of collective bands, where it 
is the coupling within the Af = 4 orbitals which induces 
the triaxial shape according to the mechanism described 
in Refs. [l3, H3|- 



A. Observed and calculated total energies 

Adding one or two neutrons, Fig. []] suggests that the 
most favoured configurations for 164 > 16 5l u w [\\ b e formed 
if these neutrons are placed in the 5/2 [523] orbital, where 
thus two bands with different signature are formed in 
164 Lu. In Fig. gfa), where the observed @, H [56| and 
calculated bands are compared, it is the lowest TSD band 
in the respective nuclei and in addition band TSD3 in 
164 Lu which are assigned to the configurations discussed 
above. Note that contrary to Ref. [56|, we have assumed 
that this TSD3 band has negative parity. The assignment 
in Ref. [H is based on Ref. [57| where band TSD3 is 
given positive parity based on the assumption that it is 
unlikely with a stretched Ml transition with such a high 
energy as 1532 keV. We find this conclusion questionable 
because in the decay of TSD1, such transitions with 1452 
keV and 1541 keV have been observed in Ref. [13] and 
Ref. (56[, respectively. Indeed, the similar decays of the 
TSD1 and TSD3 bands rather suggest to us that they 
have the same parity and this conclusion gets additional 
strong support from the comparison with calculations, 
indicating that these two bands are signature partners. 

For 163 Lu, the difference between calculations and ex- 
periment shown in the lower panel of Fig. HJa) is close 
to zero at high spin, where pairing correlations which are 
not included in the CNS formalism should be small. The 
differences are then getting larger at lower spin values, in- 
dicating the increasing importance of the pairing correla- 
tions. The curves for 165 Lu are similar leading to close to 
identical difference curves in the lower panel of Fig. [2ja) . 
The similarities between the observed bands indicate that 
the orbital which is occupied in 165 Lu but not 163 Lu is not 
strongly deformation polarizing and not giving any large 
spin contribution, as is the case for the 5/2 [523] orbital, 
which is selected in the calculations. The two bands in 
164 Lu, come close to the average of the 163 Lu and 165 Lu 
bands at high spin in Fig. HJa). Indeed, this is the case 
for all observed spin values in the (unpaired) calculations, 
while at lower spin values the odd-N energies come higher 
in experiment. This is what would be expected from a 
smaller pairing energy in the odd compared with the even 
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FIG. 2. (Color online) The observed energies of selected TSD bands for Lu isotopes are shown relative to the rotating liquid 
drop energy in the upper panels, with the calculated bands assigned to them in the middle panels and the difference between 
calculations and experiment in the lower panels. The TSD bands for 163 ' 164 > 165 l u are shown in the panels to the left and those 
of 161 ' 163 - 167 l u in the panels to the right, i.e. 163 Lu is shown in both cases to facilitate the comparison. Solid lines correspond to 
positive parity configurations and broken lines correspond to negative parity. Similarly, solid symbols correspond to signature 
a — 1/2 (a = for A even) and open symbols correspond to signature a = —1/2 (a = 1). Note how these differences are almost 
identical for most of the bands where the differences between experiment and calculations for the low spin states is understood 
from the neglect of pairing correlations in the CNS calculations. The experimental data are taken from Refs. |53| - [5S| ] 



neutron systems and it should even be possible to get an 
idea of the strength of the pairing correlations from this 
comparison. Furthermore, the calculations predict the 
correct signature for the favoured bands in 164 Lu. This 
gives additional support to the present assignments even 
though the splitting is somewhat overestimated in the 
calculations. 

Fig. Q] suggests that the additional holes in 161 Lu rela- 
tive to 163 Lu should be placed either in the «i3/2, 5/2[642] 
orbital or in the 1/2 [400] orbital. The result of the de- 
tailed calculations, see Fig.[2jb), is that the latter deexci- 
tation, i.e. the neutron configuration (42)6 is favoured for 
lower spin values while the former deexcitation, i.e. the 
neutron configuration (22)4 is favoured for higher spin 
values. Indeed, it appears that this agrees with experi- 
ment (53), [54| because in the observed band, one can see 
a smooth crossing for spin values / = 30 — 40, where the 
two unpaired configurations cross. Thus, with this as- 
signment and with our choice of spin values for the 163 Lu 
band, the difference curve in the lower panel of Fig. [2fb) 
have almost the same shape as for 161 Lu (where we have 
chosen an excitation energy of the unlinked band in 161 Lu 



similar to that for the 163 Lu band). Furthermore, with 
pairing included, the crossing between the neutron (42)6 
and (22)4 configurations will be seen as a smooth paired 
crossing within the ia.3/2 orbitals [4!|. Note that the 
two neutrons which are shifted from down-sloping to up- 
sloping orbitals lead to a considerably larger deformation 
for the (42)6 configuration, e 2 ~ 0.43, 7 ~ 23°, than for 
the (22)4 configuration, e 2 ~ 0.37, 7 - 20°. This latter 
deformation is typical for the yrast TSD bands in the 
other Lu isotopes with N = 92 — 96. 

Coming to 167 Lu, Fig. Q] indicates that the two addi- 
tional neutrons compared with 165 Lu might be put in the 
ll/2[505] orbital or in the 7/2[633] orbital. However, the 
detailed calculations show that the latter configuration 
is much less favoured for spin values above / = 30 in 
accordance with the general experience that it becomes 
energetically expensive to build spin in configurations of 
high-j shells which are half-filled or more than half-filled, 
see e.g. Fig. 12.11 of Ref. H|. As seen in Fig.^b), the 
energy vs. spin dependence of the (20)6 configuration 
in 167 Lu is close to that of the (22)6 configuration in 
165 Lu while the calculated energy is considerably higher 
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in 167 Lu than in 165 Lu in disagreement with experiment. 
This discrepancy would disappear if the ft-ii/2 subshell 
was lowered by a few hundred keV. 

There are a few more observed TSD bands in Lu nu- 
clei which we have not considered here. Thus, there are 
three unlinked bands in 162 Lu [53j]. It appears to be easy 
to assign spins and excitation energies to these bands 
so that they agree with calculations, but these assign- 
ments would be very tentative. One could note how- 
ever that the beginning of a band-crossing is observed in 
the TSD3 band which appears to be very similar to the 
band-crossing in 161 Lu suggesting a similar origin and 
thus an appreciable deformation change also in 162 Lu. 
Another band which we have not discussed here is TSD2 
in 164 Lu (56[. One could expect a neutron configuration 
with all orbitals up to the N = 94 gap in occupied but 
with a hole in the unfavoured 5/2 [642] orbital, see Fig. 
03 Indeed, the parity and signature of the observed band 
agrees with this assignment but the curvature of the E 
vs. / function of the calculated configuration appears 
too large. In addition, the observed band appears to go 
through a smooth band-crossing which is not easy to ex- 
plain. There is an interesting branch of band TSD1 at 
high spin which has a larger alignment and is referred 
to as X2 56]. This branch might be assigned to the 
configuration with the valence neutron excited from the 
favoured 5/2[523] orbital to the favoured 1/2 [770] H3/2 
orbital, see Fig. [3] below. In addition, there are several 
bands assigned as wobbling excitations in the odd Lu iso- 
topes which will of course not be described by any CNS 
configuration. 

Fig. [1] is drawn at no rotation, cj = 0, and is thus 
mainly helpful for the understanding of configurations 
at low or intermediate angular momentums. In order 
to get an understanding of the configurations which are 
favoured at higher angular momentums, it is more in- 
structive to draw a single-particle diagram at w > 0, 
which will lead to a more complicated diagram because 
the orbitals will split into two branches with signature 
a = 1/2 and a = —1/2. Such a diagram is provided 
in Fig. [3J It suggests that the favoured configurations 
for N = 90 - 94, i.e. for 161 - 165 Lu, will be about the 
same as for u = 0, but for N = 96 ( 167 Lu), it will be 
more favourable to put the two extra neutrons in the 
lowest l/2[770] orbital or in the 1/2[651] orbital (of j 15/2 
and in/2) 39/2 origin, respectively). This is also in agree- 
ment with the detailed calculations which shows that 
such a configuration becomes favoured in energy at tri- 
axial shape above I ~ 35 when combined with the same 
favoured proton configuration as for the lower spin states 
8(21). At these higher frequencies and deformations, it 
will however be favourable if also the deformation driv- 
ing second proton i 13 / 2 orbital will be occupied leading 
to the favoured 8(22) configuration for 168 Hf which will 
be discussed below. 
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FIG. 3. (Color online) Similar to Fig. [T] but with rotation 
added, i.e. a rotational frequency which increases linearly 
up to uj/ujq — 0.05 in the 'Nilsson diagram' to the left and 
then keeping this value of the rotational frequency when ax- 
ial asymmetry is added. Solid and dotted lines are used for 
positive parity and dashed and dot-dashed lines for negative 
parity, where dots are used for signature, a = — 1/2. 



B. Effective alignments, i e ff 

In our analysis of TSD bands in Lu isotopes, we will 
also consider the differences of spin, /, at a constant 
frequency, Hu), and compare the experimental and the- 
oretical data. This quantity referred to as the effective 
alignment, i e ff has been very important for the clas- 
sification of the SD bands in the A = 150 region, see 
e -g- [13) US Hil- It is a direct measure of the contribu- 
tion from different Nilsson orbitals. It is mainly useful 
when pairing can be neglected but for the Lu bands, the 
pairing correlations are rather small and we can further- 
more assume that pairing gives about the same contri- 
bution if the comparison is limited to the odd isotopes 
with an even number of neutrons. Thus, effective align- 
ments of neutron orbitals for the lowest TSD bands in 
Lu nuclei are shown as a function of rotational frequency 
(uj = E~f/2), for the experimental bands in Fig. Sta) and 
for the theoretical configurations assigned to these bands 
in Fig. |UJb). Note that in this case, i e // is a measure 
of the spin contribution from a pair of particles in the 
respective orbitals. 

The general agreement between experiment and the- 
ory in Fig. 2] indicates that we do understand which or- 
bitals are filled in the lowest TSD bands in the odd Lu 
isotopes. The spin contribution of the orbital which is 
being occupied when going from 161 Lu to the 163 Lu is 
very small and positive at Hu> < 0.5 MeV but it changes 
for fajj > 0.5 MeV where i e // turns negative. The cal- 
culated i e ff shows the same feature which can be traced 
back to a change of structure in 161 Lu from [8(21), (42)6] 
to [8(21), (22)4] at tkjj ~ 0.5 MeV. The value of i efi when 
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FIG. 4. (Color online) Effective alignment, i e ff, for Lu iso- 
topes, extracted from (a) experiment and (b) corresponding 
calculated bands. 



comparing the bands in 163 Lu and 165 Lu is close to zero 
but rather negative, corresponding to a small negative 
spin contribution from the orbital which becomes occu- 
pied. This orbital is located in the middle of the /19/2/7/2 
subshells and is labelled 5/2 [523] in Fig. Q] When two 
neutrons are added to 165 Lu, a spin contribution close to 
zero is obtained in both experiment and calculations for 
libj > 0.5 MeV. This agreement supports the assignment 
that it is the highest orbital, If /2[505], which is be- 
ing occupied. Note that this upsloping orbital will have 
a strong shape polarization, i.e. the shape ch ang e will 
have an important contribution to z e //, see e.g. [59(. The 
fact that calculations and experiment diverge at smaller 
frequencies could be caused by increasing pairing corre- 
lations so that the assumption that an orbital is either 
filled or empty is strongly violated. 

The present calculations show that the standard CNS 



formalism provides a reasonable interpretation for the 
TSD bands in Lu isotopes. However, it is questionable 
whether this approach, including approximations pointed 
out in sect. HH is suitable to study also the TSD bands 
in 168 Hf which have a larger deformation. In the next 
section, these approximations will be tested on 168 Hf 



IV. ANALYSIS OF SPECIFIC FEATURES OF 
THE CNS FORMALISM 

Representative potential energy surfaces (PES) with 
(tt, a) = (— , f ) for spins 1 = 1, 31, 41, 51, 61 are displayed 
in Fig. [5] for 168 Hf. Similar behavior is also found for the 
other (71", a) combinations. At low spins, from I = 1 to 
/ = 3f , the lowest energy minimum in the PES's corre- 
sponds to a almost prolate shape at (£2,7) ~ (0.23,0°). 
As the angular momentum increases, this minimum mi- 
grates to a somewhat larger deformation; for example 
(£2,7) ~ (0.26, 3°) at spin I = 41. For spin values 
I > 50, the minimum energy corresponds to a TSD shape 
at the deformation (e 2 ,7) ~ (0.44, 20°). 



A. Off-shell matrix elements and More shells 

As it has been pointed out in sect. [TT1 all off-shell el- 
ements in the rotating basis |n x n.2Ti3£) are small and it 
is therefore natural to neglect them. If the off-shell ma- 
trix elements are included, the shell number Af ro t will 
not the good quantum number and the rotating basis 
functions lose their advantage to diagonalize the Hamil- 
tonian matrix. It is then easier to use the stretched 
spherical harmonic basis functions \N t £tjt^t) which are 
eigenkets of the spherical harmonic oscillator Hamilto- 
nian huo{^2i l), the square of the stretched angular mo- 
mentum 2t an d its projection, jz,t- With these basis 
functions, the cranking term couples between basis state 
of the shells Aft and At ±2 which have the same signature. 

When calculating the total energy, we need the shell 
energy and the rotating liquid drop energy (Eq. (Q}). The 
addition of the off-shell elements will only effect the shell 
energy. As illustrated in Fig. [51 the shell energies ob- 
tained from the diagonalization of the Hamiltonian in the 
two cases come very close for all spin values at a large 
triaxial deformation with a typical (see below) hexade- 
capole deformation, £4 = 0.028. Note that even though 
the coupling between the Af ro t shells is neglected in the 
rotating basis, the j x term is still fully accounted for be- 
cause it is included in the basis. This is contrary to the 
stretched basis where the finite basis size corresponds to 
a (small) approximation. With more shells included, this 
approximation will be negligible. 

In the standard CNS calculations, all shells having 
the principal quantum number less than or equal to 
Af m ax = 8 are included in the diagonalization. The im- 
portant question is now if more shells are needed in order 
to reproduce the solution accurately enough for a heavy 
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FIG. 5. (Color online) Calculated potential-energy surfaces versus quadrupole deformation £2 and the triaxiality parameter 7 
of 168 Hf with (71-, a) — (— , 1) for spins 1=1, 31, 41, 51, 61. Contour lines are separated by 0.25 MeV and the 7 plane is marked 
at 15° intervals. Dark regions represent low energy with absolute minima labeled with a dot. 



nucleus like 168 Hf. Naturally, the required value of N max 
depends on particle number, the shape of the potential 
to be diagonalized and for the stretched basis also on the 
rotational frequency, uj. To illustrate the importance of 
the cut-off error, the yrast energy was calculated includ- 
ing off-shell couplings with N m ax — 12, i.e. with four 
added shells. For the specific deformation illustrated in 
Fig. [51 it turns out that the energy of the yrast line with 
Nmax — 12 does not decrease relative to the calculation 
with Nmax — 8 but it rather increases. The reason is that 
with the increase of the number of shells, both the total 
discrete and smoothed energy decrease. The total dis- 
crete single-particle energy with N m ax = 12 differs from 
that with Nmax = 8 by about 30 keV for spins I < 40 
and 120 keV for spins / > 40. Since the corresponding 
smoothed single-particle energy is shifted by about 90 
keV at spins I < 40 and 260 keV at spins J > 40, the 
resulting shell energy, 



(7) 



differs only by ~ 40 keV at spins / < 40 and ~ 140 
keV I > 40 from the corresponding value with Nmax — 
8, see Fig. [6] Thus for the equilibrium deformations of 
168 Hf in an extended spin range, the cut off at Nnax = 8 
introduces only small changes in E s h which are essentially 
negligible compared with other uncertainties. 



B. Minimization in five dimensions 

In general, for axial symmetric shapes it is only the 
£40 (with quantization around the symmetry axis) shape 
degree of freedom which is expected to be of major im- 
portance because the energy is even (independent of the 
sign) in £42 and £44. This is only valid at no rotation 
around the perpendicular axis but if the rotational fre- 
quency is not extremely high, it is still expected that only 
the £40 degree of freedom will be of major importance. 
Furthermore, shapes corresponding to small quadrupole 
deformations, are never far away from a symmetry axis in 
the (£2, 7)-plane so it should be sufficient to minimize the 
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FIG. 6. (Color online) The shell energy for the configuration 
[8(22), (22)6(11)] with deformation parameters £2 ~ 0.43,7 ~ 
20°, e 4 ~ 0.028 in the 168 Hf nucleus. The circle symbol shows 
the calculations including the off-shell elements and Nmax = 
12, the square symbol including the off-shell elements and 
Nmax = 8 and the diamond symbol excluding the off-shell 
elements and Nmax = 8. 



energy in only one £4 degree of freedom also in this case. 
This is supported by studies of the smooth terminating 
bands in 109 Sb [3, Eo] where the energy is lowered by less 
than ~ 50 keV when it is minimized in three £4 degrees of 



freedom [6l|. For a triaxial shape and large quadrupole 
deformation on the other hand, the full minimization in 
the £4i-parameter space might be more important. 

In order to make a full minimization in the five di- 
mensional deformation space, the total energy of 168 Hf is 
calculated at the following grid points: 

x = 0.18[0.02]0.44 

y = 0.08[0.02]0.42 
£ 40 = 0.005[0.01]0.045 
£ 42 = -0.02[0.01]0.02 
£44 = -0.01[0.01]0.03, 

where [x, y) are Cartesian coordinates in the (£2,7)- 
plane. The (x, y)-coordinates are connected with (£2,7) 
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FIG. 7. (Color online) The 168 Hf yrast energies relative to a rotating liquid drop energy E r jd as a function of spin / for 
the four combinations of parity and signature, (+,0), (+, 1), ( — ,0) and (— , 1). The circle symbol shows the minimum energy 
in (£2, 7, £40, £42, £44) space of deformation and the square symbol in (£2, 7, £4). Solid lines correspond to positive parity 
configurations and broken lines correspond to negative parity. Similarly, solid symbols correspond to signature a — and open 
symbols correspond to signature a — 1. The steep increase at spin values / > 72, which are most apparent for the (+,0) 
and (—,0) yrast energies is mainly caused by the rotating liquid drop reference energy which shows a discontinuity when the 
equilibrium shape moves away from the 7 = 60° axis, the superbackbend according to Ref. [3l[ (see sect. II V C|) . 



by the expressions 

x = £2 cos(7 + 30°), y = £2 sin(7 + 30°) 

In our numerical calculations, the quantization axis coin- 
cides with the rotation axis to simplify the diagonaliza- 
tion. Therefore, 7 should be replaced by (7 + 120°) in 
Eq. ©, when defining the £\i parameters. With this def- 
inition, we relabel the principal axis but the same nuclear 
shapes are formed in the £4i-space. Especially, it is for 
rotation around the symmetry axis ((7 = 60, —120°) that 
axially symmetric shapes are formed with only £40 7^ 
while axially symmetric shapes at 7 = 0° are described 
by all e 4 i 7^ 0. 

In Fig. the 168 Hf yrast energies are drawn relative 
to a rotating liquid drop energy E r ;^ as a function of 
spin / for the four combinations of parity and signa- 
ture, {it, a) = (+,0), (+, 1), (-,0) and (-,1). They are 
compared with the corresponding energies from the min- 
imization in the (£2, 7, £4) parameter space. In our cal- 
culations, the reference energy E r id is minimized in a 
deformation space (£2, 7, £4) for each spin value. 

As one can see, at spins 10 < J < 45, the yrast states 
in the deformation space (£2, 7, £40, £42, £44) are only a 



few keV lower in energy than that of in the (£2, 7, £4) de- 
formation space. On the other hand, the gain in energy 
in the high spin region, / > 45, is important and amounts 
to 0.5 MeV at some spin values. These findings are con- 
sistent with the general expectations discussed above. 
Thus, according to the potential-energy surfaces in the 
CNS calculations for 168 Hf (see Fig. [5]), the yrast states 
are built from configurations which have prolate shape 
with £2 ~ (0.23 — 0.26) for spin values below / ~ 45 
but at non-axial shape with (£2,7) ~ (0.44,20°) (TSD 
shapes) for spins I > 45. Therefore in the following, 
we do the minimization process in the deformation space 
(£2, 7, £4) to study the bands close to axial shape and 
in the deformation space (£2, 7, £40, £42, £44) to study the 
TSD bands in 168 Hf. 

In order to illustrate the variation of the £4^ parameters 
in the two cases, they are drawn in Fig. [S] as functions 
of spin I for the TSD configuration, [8(22), (22)6(11)]. In 
the complete minimization, the £4^ parameters get dif- 
ferent values relative to Eq. ([6]) in the full spin range, 
/ = 30 — 80. The value of the £40 parameter becomes con- 
siderably larger, £40 ~ 0.035 compared with £40 ~ 0.020 
in the restricted variation. The £42 parameter changes 
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sign over most of the spin range while the £44 parameter 
varies faster and gets larger values. 

The discontinuity in the variations of the en param- 
eters at spin / ~ 65 is understood from a crossing of 
high-j and low-j orbitals in this configuration which is 
explained below. The energy surfaces at spin I = 50 and 
for the same [8(22), (22)6(11)] configuration are shown 
in Fig. |9fa-c), in the planes (£40,242), (£42, £44) and 
(£44, £40) for a constant value close to the minimum of 
the third parameter. These figures indicate that the to- 
tal energy is well-behaved with only one minimum in the 
(£40, £42, £44) space. 



C. The reference energy 

In order to highlight the details of high-spin bands, 
their energy is often shown relative to a reference. For 
a long time, the standard choice of such a reference has 
been E re f = C ■ 1(1 + 1) MeV/ti 2 , where C is a constant 
23] for a specific nucleus. In calculations based on the 
CNS approach, the constant has generally been chosen 
as C = 32.32A 5 / 3 MeV [H], which means that the ref- 
erence energy corresponds to rigid rotation at a prolate 
deformation, e = 0.23, assuming a sharp nuclear radius 
tqA 1 / 3 with tq = 1.2 fin. With this choice, the increase 
or decrease of E(I) — E re f is relevant and it becomes in- 
structive to compare rotational bands in different mass 
regions. On the other hand, different constants have been 
used in the literature so one should be careful before 
drawing any conclusions from the slope of E(I) — E re f 
curves. For examples, while the A-dependent expression 
specified above gives C = 0.00665 for A — 163, the value 
C = 0.0075 has often been used for the TSD bands in Lu 
nuclei, see e.g. Refs. |53 - l56| . This larger value of C leads 
to a substanaial down-slopes for the observed energies of 
these bands, while these energies are rather constant with 
our standard choice for C. 

The absolute value of E(I) — E re f is dependent not 
only on the [shell] energy for a specific spin value but also 
on the [shell] energy at the ground state. This appears 
reasonable for low- and intermediate-spin states formed 
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FIG. 9. Energy surfaces shown as functions of two of the 
three £4; parameters for the configuration [8(22), (22)6(11)] at 
spin I = 50 and quadmpole deformation parameters (£2,7) = 
(0.43, 20°). The contour line separation is 0.1 MeV in (a) and 
0.2 MeV in (b) and (c). 



at similar deformation as the ground state. However, for 
higher spin values, the deformation or coupling scheme 
can be quite different and it is then more reasonable to 
find an absolute reference, independent of the ground 
state for that specific nucleus. Such an absolute reference 
is provided by the rotating liquid drop (RLD) model [62j , 
which can be used in a similar way as a static liquid drop 
model is used for nuclear ground states [HI, With 
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this in mind, a RLD reference was introduced in Ref. [25| , 
where it was concluded that a good fit to nuclear high- 
spin states could be achieved using the Lublin-Strasbourg 
drop (LSD) model for the static liquid drop energy 
with the rigid body moment of inertia calculated with 
a radius parameter ro = 1.16 fm and a diffuseness a — 
0.6 fm (63. With this choice, it becomes possible to 
describe the absolute energy of nuclear high-spin state 
with a similar accuracy (~ ±1 MeV) as nuclear masses 
251. 
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FIG. 10. (a) The difference between the reference energy 
based on the LSD model with a moment of inertia calculated 
from a diffuse surface [2^], and the standard 1(1+1) reference 
for the 168 Hf nucleus/ (b) The 168 Hf yrast energies relative 
to standard 1(1 + 1) reference. These energies, which are 
minimized in the (£2,7,24) space of deformation, cf. Fig. [7] 
are shown as a function of spin / for the four combinations of 
parity and signature, (+,0), (+, 1), ( — ,0) and ( — , 1). 

The rotating liquid drop energy at its equilibrium 
deformation is plotted relative to the fixed reference 
C ■ 1(1 + 1) in Fig. [LOTah This value is thus showing 
the difference what concerns spin dependence of the 'pre- 
vious' and 'present' reference energies. Note that both 
these references are the same for all bands in one nu- 
cleus, but that the mass dependence is somewhat differ- 
ent. It is easy to understand the general structure of the 
curve in Fig. 110( a). At low spin values, the equilibrium 
deformation of the rotating liquid drop energy is spheri- 
cal corresponding to a small moment of inertia and thus 



a larger reference energy. With increasing spin, the in- 
creasing oblate deformation of the rotating liquid drop 
energy corresponds to an increasing rigid body moment 
of inertia and at I sa 50, the difference starts to decrease 
corresponding to the same moment of inertia for the two 
reference energies. At even higher spin at I « 74, the 
so-called superbackbend occurs [3l], |6(| , when the rotat- 
ing liquid drop energy loses its stability towards triaxial 
shape. This corresponds to a rapid increase of the rigid 
body moment of inertia, leading to large negative values 
for higher spin values in Fig. HUT a) . 

It is now easy to understand the differences when the 
yrast energies are plotted relative to the two differences in 
Figs. 171 and [POlbh respectively. Thus the general appear- 
ance is the same up to I — 60 — 70 but with a larger ten- 
dency for decreasing values at low spin with the rotating 
liquid drop reference. The large differences are however 
at the highest spin values where the equilibrium defor- 
mations in the CNS calculations are generally found at a 
large deformation with a small moment of inertia which 
corresponds a large down-slope when this energy is shown 
relative to the C -1(1+ 1) reference, see Fig.fTOTb). With 
the rotating liquid drop reference on the other hand, the 
reference energies and CNS energies will on the average 
have the same spin dependence but a not so nice feature 
is that the large changes in the reference energy at the 
superbackbend leads to a somewhat strange behaviour of 
the energies at 7 « 74 in Fig. [7] 

Let us also point out that the smaller radius parameter 
combined with the diffuseness correction corresponds to 
essentially the same rigid moments of inertia in the two 
reference energies for mass numbers A = 150 — 200. For 
smaller mass numbers on the other hand, the diffuseness 
correction becomes more important. For example, in the 
A = 60 region, the spin dependence of the two references 
is very similar for spin values I = — 15 but they become 
quite different at higher spin values. Thus, already at 
I = 30, the energy of the rotating liquid drop reference 
is 2-3 MeV smaller than the standard C' 1(1+1) reference. 



V. THE HIGH-SPIN BANDS IN 168 HF 

A. Observed high-spin bands in 168 Hf 

Experimental excitation energies relative to a rotat- 
ing liquid drop energy, E r ;d, as a function of spin I and 
spin, kinematic and dynamic (J^) moment of in- 

ertia as a function of rotational frequency, hw, are drawn 
in Fig. UTTa-d). respectively, for the five bands in 168 Hf 
which are observed well beyond I = 40, where pairing 
correlations should be negligible. From Fig. [TTl a) one 
can see that there is a break in the rotational pattern at 
I ~ 12 and 7-40 in band 1 and at I ~ 20 and 7-40 in 
band 3. Furthermore, the spin 

(Fig. HEM) and the J« 
moment of inertia (Fig. Illf c)) are triple- valued for band 
1 at I ~ 12, i.e. band 1 goes through a full backbend 
at this spin value. The source of this backbend is the 
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FIG. 11. (Color online) (a) Experimental excitation energies relative to that of a rotating liquid drop E r ;d as a function of 
spin J, (b) spin as a function of the rotational frequency, (c) kinematic (J 1 - 1 ') and (d) dynamic (J' 2 ') moments of inertia as 
a function of rotational frequency, tiw for band 1, band 3, ED, TSD1 and TSD2 bands in 168 Hf. The data are taken from 
Refs. UMl. 



decoupling and spin alignment of an i\3/ 2 neutron pair 
from the pairing field [671 ]. The unsmoothness in and 
a small peak in (Fig. UHd)) at I ~ 20 in band 3 in- 
dicates a weak crossing at this spin. The larger variation 
of and a huge jump in at I ~ 40 (tvjj ~ 0.55) 
correspond to a larger spin alignment in band 1 and band 
3 at this spin. The excitation energy varies smoothly for 
the ED, TSD1 and TSD2 bands, which means there is 
no crossing in these bands, even though the ED band 
displays a small rise or bump in the value with the 
maximum at ~ 0.45 MeV. 



B. Calculated rotational band structures in Hf 

For the prolate shape minimum (at e 2 ~ 0.23) for 
(•7T, a) — (+,0) and (— , 1), calculated excitation energies 
for the low-energy configurations in 168 Hf are plotted, rel- 
ative to that of a rotating liquid drop in Figs. H2T a) and 
(b) , respectively. As pointed out in section IIVB1 these 
configurations are obtained from energy minimization in 
the deformation space (£2, 7, £4)- The (+,0) yrast line 
has the configuration 7r(/i 11 / 2 ) 8 ^(«i3/2) 4 or [8,4] in the 
shorthand notation for spins 1^0 — 40, see Fig. H5T a). 
As the angular momentum increases, the lowest state 
is obtained by exciting a proton from a high-j orbital 
of hx\/2 character to an orbital of /19/2/7/2 character in 
the Af = 5 shell. Therefore the yrast line is built from 



the Tr(h 11/2 ) 7 {h 9/2 f 7/2 ) 1 v(i 13/2 ) i orbitals or [7(10),4] in 
a short spin range for / > 40. Then for / ~ 44, the calcu- 
lated yrast configuration is [8(11), 5] before the [7(11), 4] 
configuration comes lowest in energy at I ~ 50. The 
single-particle occupancy in these configurations can be 
understood from Figs. ITST a) and fTBT b). where the single- 
particle Routhians are plotted for protons and neutrons, 
respectively. The configuration change in the (+, 0) yrast 
states at I ~ 40 is explained from the crossing between 
the 7/2 [523] and 1/2 [541] orbitals at hw ~ 0.6 MeV. 
There is a large single-particle shell gap associated with 
neutron number N = 96 that continues to Huj ~ 0.7 MeV, 
see Fig. [TBTb). so the neutron configuration v(ii3/ 2 ) 4 is 
favoured up to spin values beyond / = 50. 

The study of the calculated excitation en- 
ergies for the low-energy configurations with 
(•7T, a) = (— , 1) and axially symmetric shapes 
(Fig. [T2Tb)) suggests that the (—,1) yrast line is 
built on configurations [8(10), 4] and [8(11), 4] which 
correspond to 7r(/i 11/ 2) 8 (/i 9 /2/ 7 /2) 1 ^(*i3/2) 4 and 
7i"(/in/2) 8 (^9/2/7/2) 1 (*i3/2) 1 ^(«i3/2) 4 , respectively. 
These two bands which cross at / ~ 40 have the 
same («i3/2) 4 neutron configuration. The calculated 
deformations are (e 2 ,7) - (0.23,0°) and (0.26,5°) 
for [8(10), 4] and [8(11), 4] configurations, respectively. 
In Fig. [CUb), also the [8,5] and [7(10),5] configu- 
rations are drawn. They have normal deformation, 
(£2,7) ~ (0.23,0°) and cross at spin / ~ 50. The config- 
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FIG. 12. (Color online) Calculated total excitation energies 
in 168 Hf relative to a rotating liquid drop reference for low- 
energy configurations with (a) (it, a) = (+, 0) and 7 ~ 0°, (b) 
(it, a) — (— , 1) and 7 ~ 0°. Each band is shown by a label 
which is explained in text. 



uration 7r(ft 11/2 ) 8 (/i9/2/7/2) 1 («i3/2) 1 ^(^ii/2) 2 (*i3/2) 6 or 
[8(11), (02)6] in shorthand notation has the deformation 
(£2,7) ~ (0.3,1°). In fact, the two holes in /in/2 
neutron orbitals (see Fig. fT3Tb)) lead to an enhanced 
deformation. 

The calculated energies at the TSD minimum are 
drawn in Fig. 114( a) for six low-energy configurations of 
168 Hf. The associated dynamic moments of inertia are 
given as a function of rotational frequency in Fig. 114T b) . 
As discussed above, the total energy is minimized in a 
five dimensional deformation space (£2, 7, £40, £42^44) in 
this case. All TSD bands are built on the proton con- 
figuration n(h 11/2 ) s (h 9/2 f r/2 ) 2 (i }3/2 ) 2 or [8(22)]. This 
is understood from a proton single-particle shell gap 
at (£2,7) ~ (0.43,20°) for Z = 72 which is seen in 
Fig. rrSF a). For neutrons at TSD deformation, a large 
energy gap is calculated for N — 97 as anticipated from 
Fig. H and seen Fig. H5Tb). This suggests that 169 Hf 
should be a good candidate to observe TSD bands ex- 
perimentally. The lowest N = 96 configurations are 
formed from a neutron hole in the two signatures of the 
(57/2^5/2)5 (^9/2/7/2) and (i u /259/2) orbitals below the 
N = 97 gap, resulting in six different neutron configu- 
rations. Five of these together with one configuration 
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FIG. 13. (Color online) Single-particle proton (a) and neutron 
(b) energies as a function of rotational frequency (Routhians) 
at the deformation £2 ~ 0.23, 7 ~ 0° and £4 = 0.009. The 
orbitals are labeled at hw = by the asymptotic quantum 
numbers. A few important orbitals for the present inter- 
pretation are labeled by their dominating j-shell. The line 
types distinguish between different (tt, a) combinations: solid 
lines represent (+, +1/2), dotted lines (+, —1/2), dashed lines 
(-, +1/2) and dash-dotted lines (-, -1/2). 



with two Af = 7 neutrons are combined with the favored 
proton configuration, forming six the low-energy triaxial 
structures shown in Fig. H"4T a). Note that all the neu- 
tron configurations are built on six neutrons in 113/2 or- 
bitals, two holes in hu/2 orbitals and two, three or four 
holes in AT = 4 orbitals. The calculated deformation is 
(£2,7) ~ (0.49,22°) for the [8(22), (42)6(22)] configura- 
tion while it is (£2,7) ~ (0.4 - 0.45,20°) for the other 
TSD configurations. All of the theoretical TSD bands 
shown in Fig. IT47 a) display a decreasing value of J*- 2 ) 
with increasing rotational frequency. At hw < 0.45 MeV, 
the value of J 1 - 2 ' decreases more strongly for the con- 
figuration [8(22), (42)6(22)] and more smoothly for the 
[8(22), (22)6(11)]. The values of are very close to- 
gether at hu > 0.45 MeV and only the configuration 
[8(22), (22)6(11)] experiences a sharp discontinuity in the 
j( 2 ) moment of inertia at hw ~ 0.7 - 0.8 MeV. This 
discontinuity is because of a crossing in the neutron sin- 
gle particle orbitals between a high-j 213/2 orbital and 
a low-j (in/259/2) orbital at hw ~ 0.7 MeV. The cross- 
ing is indicated by a circle in Fig. fTBTb). In the other 
TSD configurations the (in/259/2) orbital has been filled 
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FIG. 14. (Color online) (a) Calculated total excitation ener- 
gies relative to a rotating liquid drop reference as a function 
of spin for six low-energy configurations with TSD shape and 
(b) J (2) values as a function of rotational frequency for the 
six low-lying collective configurations in 168 Hf. Each band 
is shown by a label which is explained in text. Solid lines 
correspond to positive parity configurations and broken lines 
correspond to negative parity. Similarly, solid symbols cor- 
respond to signature a — and open symbols correspond to 
signature a = 1. 
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FIG. 15. (Color online) Single-particle proton (a) and neutron 
(b) energies as a function of rotational frequency (Routhians) 
at the deformation £2 ~ 0.43, 7 ~ 20° and £4 = 0.03. The 
orbitals are labeled at huj — by the TV shell to which they 
belong with the ordering within the M shell as a subscript. A 
few important orbitals for the present interpretation are also 
labeled by their dominating j-shell(s). The line types distin- 
guish between different (71", a) combinations: solid lines repre- 
sent (+,+1/2), dotted lines (+,-1/2), dashed lines (-,+1/2) 
and dash-dotted lines ( — , —1/2). 



Comparison between calculated and 
experimental bands in 168 Hf 



and therefore the strong alignment at Hlu ~ 0.7 MeV is 
blocked and there is no anomaly in the moment of 
inertia of them. 



In the calculations, no distinction is made between low- 
j and high-j orbitals at this large deformation, i.e. only 
the number of particles of signature a = 1/2 and a = 
— 1/2 in each Av t-shell is fixed. On the other hand, 
the configurations are labeled as if such a distinction is 
made. The labels in Fig. H~4T a) refer to the configuration 
for spin values below I ~ 60. For example, the energy 
of the configuration labeled [8(22), (22)6(11)] comes down 
at spin I ~ 62, because of the crossing discussed above. 
Thus, the band should be labeled [8(22), (22)5(21)] for 
higher spin values. 



In the upper panels of Fig. 1161 experimental excitation 
energies relative to a rotating liquid drop energy for band 
1, band 3, the ED band and the TSD1 and TSD2 bands 
are drawn as a function of spin. The middle panels of 
Fig. [16] displays the calculated bands which seem to be 
closest to these experimental bands. In the lower pan- 
els, experimental and theoretical bands are compared (in 
attention to their parity and signature) and their differ- 
ences are illustrated. 



1. Band 1 

The observed band 1 has a positive parity and signa- 
ture a = 0. Therefore one can find out the structure 
of this band from the search among the lowest-energy 
configurations which have (tt, a) = (+,0) (Fig. HUT a)). 
As it is pointed in section IV Bl in Fig. HUJa) the [8,4] 
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FIG. 16. (Color online) Experimental energies (top panels) 
and theoretical energies (middle panels) relative to a rotating 
liquid drop and their differences (lower panels) as a function 
of spin for (a) band 1 and band 3 and (b) ED, TSD1 and 
TSD2 bands in 168 Hf. 



and [7(10), 4] configurations are the lowest states in en- 
ergy at spins I < 40 and I > 40, respectively. The [8,4] 
configuration has an even number of neutrons in 
orbitals. Thus the observed backbending at / ~ 12 in 
band 1 (see section IV A[) could occur in this configura- 
tion. Furthermore, the change of structure from [8,4] to 
[7(10) ,4] which happens at spin / ~ 40 corresponds to 
the observed break in the rotational pattern at / ~ 40 in 
this band (see Figs. [TlTa-d)). Band 1 is compared with 
the [8,4] and [7(10), 4] configurations in the lower panel 
of Fig, [TBTa). As one can see the differences between the 
theoretical and experimental data for band 1 is rather 
constant and about —1 MeV at spins / > 35 if a transi- 
tion occurs from [8,4] to [7(10),4]. 

The same structure has been obtained for band 1 in 
Ref. [HI using the Ultimate Cranker code [29], [68| . How- 
ever, it is concluded [l6| that the occupation of the 1/2 
[541] orbital at I ~ 40 is related to the crossing be- 
tween the proton orbitals 9/2[514] and 1/2 [541]. This 
is contrary to our calculations, see Fig. ITBTa). where the 
9/ 2 [5 14] orbital is above the Fermi surface and thus the 
transition is because of the crossing between the 7/2 [523] 
and 1/2[541] orbitals. A closer look at Fig. 5 of Ref. pj| 
indicates that there are two hu/2 quasiparticles at simi- 
lar energies where one should then mainly correspond to 
a hole in the 7/2 [523] orbital and the other to a particle in 
the 9/2 [514] orbital. Then, it appears that if the number 
of particles should not be changed drastically, the added 
particle in 1/2 [541] should be combined with a hole in 
7/2[523], contrary to the conclusion in Ref. [161 ]. 



2. Band 3 

Band 3 is observed from I = 7 to I = 47 and has 
(-7T, a) = (—,1). Excitation energies, spin, jW and 
behavior of this band are close to those of band 1 (see 
Figs. [TlTa-d)). Thus it seems that this band has a simi- 
lar deformation as band 1. The lowest-energy states with 
(7r,a) = (-, 1), see Fig.Q^b), suggest that the [8(10),4] 
and [8(11), 4] configurations should be assigned to band 
3. However, as we see in Figs. ITTT a-c). in contrast to 
band 1, band 3 does not backbend at low spins. There- 
fore the neutron configuration of this band could not be 
the same as that of band 1 (which has four neutrons in 
H3/2)- As pointed in section IVBl the [8,5] and [7(10), 5] 
configurations have almost the same deformation as band 
1 and they also have an odd number of neutrons in H3/2 
orbitals. Thus, even though these two configurations are 
calculated about 0.5 MeV higher in energy than [8(10), 4] 
and [8(11), 4], see Fig. rT2T b). they are more suitable can- 
didates for band 3. As one can see in Fig.[T27b). there is a 
crossing between [8,5] and [7(10), 5] at spin I ~ 50 which 
is in agreement with the observed crossing in band 3 ex- 
perimentally. In the lower panel of Fig. [T67a) band 3 has 
been compared to the [8,5] and [7(10), 5] configurations. 
With the transition from the [8,5] to the [7(10), 5] config- 
uration, the differences between the calculated bands and 
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band 3 is almost constant at -0.5 MeV for spins I > 35. 
Therefore band 3 is built from the («i3/2) 5 neutron config- 
uration and the proton configuration is the same as that 
for band 1, (h ll/2 f at / < 50 and (h 11/2 )' 7 \h 9 / 2 f 7/2 ) 1 for 
I > 50. This interpretation is similar to that of Ref. jl6| . 
but with the same difference as discussed above for band 
1. 



3. Band ED 

The ED band (called TSD2 in Ref. [3) has (n,a) = 
(— ,1) and is observed from I = 23 to J = 49. 
The calculations, as depicted in Fig. [T2Tb). show that 
^(/iii/2) 8 (^9/2/7/2) 1 (H3/ 2 ) 1 ^(«13/2) 4 or the [8(11),4] con- 
figuration is lowest in energy for spins / > 40 This config- 
uration, which corresponds to axially symmetric shape, 
has been suggested for the ED band in 168 Hf 15]. The 
calculated quadrupole deformation value, £2 ~ 0.26,7 ~ 
5°, is near normal deformation and far from that of 
the ED bands in the other Hf isotopes, £2 ~ 0.3 jicj . 
The suggested configurations for the ED band in Hf iso- 
topes are all built on the same proton configuration, 
*i3/ 2 ^9/2i but they are coupled to different neutron con- 
figurations jig. Ha]. 

A common feature of most interpretations of strongly 
collective bands is that only the high-j intruder orbitals 
from the higher shells are listed explicitly in the config- 
urations. These high-j orbitals are important to build 
the spin but on the other hand, it is rather the extruder 
orbitals from the lower shells which build the collectiv- 
ity. This is evident for the smooth terminating bands [24| 
and it has been underlined that it is the case for more 
collective bands, see e.g. [13, HI!- However, to our knowl- 
edge, it is only in the present CNS approach that meth- 
ods have been developed to fix the number of particles 
in the extruder orbitals. In the present case, the high- 
est hn/2 orbital is just below the Fermi surface. Thus 
we consider the configuration with two holes in the up- 
sloping ll/2[505] orbital, i.e. [8(11), (02)6]. This con- 
figuration is about 1 MeV higher than [8(11), 4] in en- 
ergy, see Fig. IT2Tb). However, as pointed in section IV B[ 
the [8(11), (02)6] configuration has a larger value for the 
quadrupole deformation, (£2,7) ~ (0-3,1°), which is in 
agreement with that of in the other Hf isotopes. Espe- 
cially, the experimental properties of the ED band are 
clearly different from those of the valence space band 
and it is only with the holes in the neutron orbitals 
that also the theoretical configuration becomes clearly 
different from the valence space configurations. With 
the [8(11), (02)6] interpretation for the ED band, the 
difference between calculations and experiment becomes 
small and almost constant as seen in the lower panel of 
Fig.lTBTb). It also suggests that such configurations with 
holes on the /111/2 neutron orbital, below N = 82, should 
be investigated for the ED bands in other Hf isotopes. 
Furthermore, as mentioned in Ref. jl7j . it appears that 
the same mechanism with holes in the hn/2 neutron or- 



bitals is responsible for the large quadrupole moment in 
the SD band of 175 Hf. 

One problem with the present interpretation is that the 
[8(11), (02)6] configuration is calculated at an excitation 
energy which is somewhat too high relative to the config- 
urations assigned to band 3. However, these differences 
are clearly within the expected uncertainties. For exam- 
ple, if the neutron h\\n subshell was placed 0.5 MeV 
higher in energy, the [8(11), (02)6] configuration would 
be calculated close to yrast. Note also that when the 
[8(11), (02)6] configuration is not calculated as yrast, it 
is straightforward to study it only in approaches like the 
present one where it is possible to fix the number of holes 
(or particles) in specific orbitals. 



4. TSD1 and TSD2 Bands 

The observed TSD bands in 168 Hf [lj| have not been 
linked to the normal-deformed level scheme so their spin 
and excitation energy can only be estimated. The struc- 
ture 7r(«i3/ 2 ) 2 ^(ji5/2«i3/ 2 ) has been suggested as the 
most probable intrinsic configuration for the TSD1 band. 
Our calculations with a complete minimization in £4 show 
that the [8(22), (22)6(11)] configuration is the lowest- 
energy TSD configuration at spins / < 50, see Fig. H4T a') . 
Note that these two suggested configurations are identi- 
cal with the same high-j orbitals occupied but, in the 
unpaired CNS formalism, also the occupation of other 
orbitals are specified, including the extruder neutron or- 
bitals with their main amplitudes in the /iix/2 and J\f = 4 
shells, respectively. 

The [8(22), (22)6(11)] configuration with (tt, a) = 
(— ,0) is about 0.5 MeV lower than the next lowest- 
energy TSD configuration ([8(22), (32)6(21)] with a = 1) 
at spins I < 45. Band TSD1 has been measured up to 
7-60 where the [8(22), (22)6(11)] configuration is cal- 
culated only a few hundred keV above the yrast line. 
Thus, we choose this configuration as our favored candi- 
date for TSD1 and compare the two bands in lower panel 
of Fig.fTBTb'). This configuration choice suggests that the 
observed band TSD1 has negative parity and even spin 
(a = 0). 

TSD1 is plotted with an assumed bandhead spin 
of I = 34h and with an energy of 11.6 MeV for 
the bandhead. This leads to a good fit between 
the observed band and the calculated configuration 
[8(22), (22)6(11)], where the shorthand notation cor- 
responds to n(h 11/2 ) s (h 9/2 f 7/2 ) 2 (i 13/2 ) 2 and v(Af = 

4) _2 (/iii/ 2 ) _2 (ii3/2) 6 («ii/259/2) 1 (ji5/2) 1 for the occupa- 
tion of the open proton and neutron sub-shells. The 
calculated quadrupole moment of the configuration 
[8(22), (22)6(11)] is in the range 11.8-9.9 eb for spin values 
/ = 30 — 80 which is in agreement with the experimen- 
tally measured value of Q t = 11.4l*;2 eb [13]. However, 
the calculated quadrupole moment is similar for all con- 
figurations in the TSD minimum so it does not help to 
discriminate between the different configurations listed 
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FIG. 17. (Color online) Theoretical and experimental dy- 
namic moments of inertia j' 2 ^ as a function of rotation fre- 
quency huj for 168 Hf. 



above. 

The configuration assignment to the TSD1 band is 
certainly preliminary and one might for example ar- 
gue that we should rather choose a configuration which 
is calculated yrast at the highest observed spin value, 
I - 60. This would then rather suggest [8(22), (22)6(21)] 
or [8(22), (32)6(21)] as the favoured choice, i.e. configu- 
rations with one more neutron excited to the in/299/2 
orbitals. In addition to the high-j particles, the TSD 
minimum is characterized by at least 2 hn/2 and 2 Af = 4 
neutron holes. This is the case also for the calculations 
presented in Ref. [l5[ even though these holes are not 
specified in the configuration labels used in that refer- 
ence. As pointed out in section Hill these holes are im- 
portant to create the smooth collective bands where it is 
mainly the Af = 4 holes which induces the triaxial shape. 

Based on comparisons with calculations and with the 
TSD1 band, see Figs. IT47 a) and [TfjT b). the spin and 
bandhead energy for the TSD2 band are estimated to be 
41H and 14.7 MeV, respectively. As pointed out above, 
the next lowest TSD configuration for I ~ 40 — 50 is 
[8(22), (32)6(21)] with (tt, a) = (-,1). This configura- 
tion is yrast for 55 < I < 70. Therefore it seems 
that this configuration is a reasonable candidate for the 
observed TSD2 band. As one can see in the lower 
panel of Fig. fTBT b). the energy difference between the 
TSD2 band and the [8(22), (32)6(21)] configuration is 
small and rather constant for I > 50. Considering the 
configurations in Fig. IT47 a). another possible choice is 
the [8(22), (22)6(21)] configuration with (n,a) = (+,1). 
Thus, the present calculations suggest that compared 
with TSD1, TSD2 has the same proton configuration 
but with one neutron excited to «n/2<?9/2 from either the 
Af = 4 orbitals or from the /19/2/7/2 orbitals. This leads 
to odd spin values (a = 1) but undetermined parity for 
TSD2. 

The configurations of the TSD1 and TSD2 bands could 
be interpreted by considering the behavior of the dy- 



namic moment of inertia. Although the spin assignments 
for these TSD bands may need revising, the dynamic mo- 
ments of inertia are not affected by these changes. A 
smooth decrease in the J*- 2 - 1 moment is observed for the 
TSD1 and TSD2 bands, see Fig. [TTl which is consistent 
with t he g eneral trend of TSD bands in the other Hf iso- 
topes UHil, Iz3]- Fig. [i7] also displays three configura- 
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tions that have characteristics similar to the two observed 
TSD bands. On the other hand, the absolute value of 
is somewhat smaller in calculations than in experiment, 
which can also be concluded from the positive curvature 
in the difference curves in the lower panel of Fig. |T5fb) . 
The TSD1 band and the [8(22), (22)6(11)] configuration 
have rather similar slopes in throughout the observed 
frequency range. The value of for the TSD2 band 
has a behavior similar to that of two suggested configura- 
tions but the calculated dynamic moment of inertia is the 
same for the others TSD configurations, see Fig. [T4Tb). 
so it does not help to choose a favorable configuration for 
the TSD2 band. 



VI. ADDITIONAL COMMENTS ON THE TSD 
BANDS IN A = 158 - 168 NUCLEI. 

A. Full minimization in Lu isotopes and 158 Er 

We have examined the full minimization approach in 
the hexadecapole deformation space (see sect. IIVB[) for 
TSD configurations in Lu isotopes. Our calculations 
show that this effect will typically decrease the minimum 
energy by 200 keV in the observed spin region (7 < 50). 
The maximum gain of about 300 keV is obtained in 
165 Lu, which has a large hexadecapole equilibrium defor- 
mation, £4 ~ 0.04. These effects will lead to some minor 
corrections on the results presented in sect. Mil but they 
will clearly not change the general conclusions. 

Similar calculations have been carried out for 158 Er 
where, according to studies in Ref. [2(1 HJ, three 
well-defined TSD minima with deformations (£2,7) ~ 
(0.37, ±20°) and (£2,7) ~ (0.45,20°) are seen. Our cal- 
culations show that a complete minimization in £4 has 
only a small influence on the energy. The gain in energy 
is always smaller than 200 keV for the TSD configura- 
tions in 158 Er. 



B. Effective alignments in a larger mass range from 



"Er to 



'Hf 



It is instructive to consider the alignment in a larger 
mass range outside the Lu isotopes. Therefore, in Fig. 1181 
the difference in alignment between the lowest TSD 
bands in 163 Lu and 168 Hf is shown for the experimen- 
tal bands and for the configurations we have assigned 
to these bands. In the calculations, we have chosen the 
spin values for the unlinked band in 168 Hf based on the 
present calculations and this will naturally lead to a gen- 
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eral agreement between 163 Lu and 168 Hf what concerns 
the effective alignment. Considering the orbitals which 
are filled in 168 Hf but not in 163 Lu, it is mainly the H3/2 
proton and the j'15/2 and in/2 neutrons which build the 
large effective alignment of almost 12 h while the filling 
of two M = 5 neutron orbitals will only have a small 
contribution to i e ff- 

It is then also instructive to compare the spin difference 
between the lowest TSD bands in 158 Er and 163 Lu which 
is drawn for experiment and calculations in the upper 
and lower panels of Fig. [TBI respectively. Several possible 
theoretical assignments are shown, corresponding to the 
lowest-energy configurations in the minima with £2 ~ 



0.34,7 « 20° (TSD1), £ 2 « 0.34, 7 w -20° (TSD2) and 
e 2 ~ 0.43,7 fa 25° (TSD3), where the different minima 
are labelled as in Ref. [21|. Furthermore, it should be 
noted that the spin values in 158 Er are not known and 
have been chosen in the range I = 23— 65ft as suggested in 
Ref. [TtJ ■ If these spin values are increased (decreased) by 
AI ft, it will correspond to a constant decrease (increase) 
for values of the curve in the upper panel by AI ft, but 
with no change of the spin dependence. 

The TSD1 configuration of 158 Er and the TSD band 
of 163 Lu have similar deformations so the corresponding 
value of i e f / measures the spin contribution from the or- 
bitals which are filled in 163 Lu but not in 158 Er, i.e. 2 
/in/2 and one hg/2 protons and 2 113/2 neutrons. They 
will then give a negative contribution to i e ff in agree- 
ment with general expectations for orbitals in the middle 
of a j-shell, see e.g. Fig. 27 of Ref. (23|. When it comes 
to the configurations in the TSD2 and TSD3 minima, the 
value of i e ff does not correspond to the contribution of 
any specific orbitals because these configurations do not 
have a common core with 163 Lu. 

In any case, it is still possible to define the difference in 
spin value for a fixed frequency and the comparison be- 
tween experiment and calculations in Fig. 1181 shows that 
it is necessary to increase the spin values in 158 Er by 
4 — 8ft to get agreement at the highest frequencies for the 
TSD3 configurations. For the TSD2 configuration on the 
other hand, experiment and calculations come close for 
all frequencies with present spin values. Note especially 
that a down-slope is seen both in experiment and calcu- 
lations for frequencies fku w 0.35 — 0.50 MeV. This down- 
slope corresponds to an additional alignment of ~ 4ft in 
this frequency range for 158 Er relative to 163 Lu. Such 
an alignment gives rise to a bump in the J*- 2 -* moment 
of inertia as discussed in some detail for the correspond- 
ing bands in 159 . 160 Er in Ref. [2(3]. As discussed there, 
the alignment is caused by a crossing between an ftn/2 
proton orbital and the lowest hg/2 orbital. As seen in 
Fig. 1181 while the observed alignment is approximately 
reproduced for the TSD2 configuration, no similar align- 
ment is seen in the TSD1 and TSD3 configurations. In 
Ref. [2(| , a bump in the moment of inertia for the 
TSD1 band, caused by a crossing between the Af = 6 
neutron orbitals, was discussed. This is however a con- 
siderably broader bump corresponding to an alignment 
in a larger frequency range which is not seen as any well- 
defined alignment in Fig.[l8j For the TSD3 configuration, 
no specific alignment is seen in Fig. [TS] and no crossing 
between orbitals is observed in Fig. [T51 which could give 
raise to such an alignment. 

The conclusion from the present analysis of the align- 
ments would thus be that the TSD2 configuration, i.e. 
the 7 < configuration, which is the preferred assign- 
ment for the yrast TSD band in 158 Er. This is also the 
preferred configuration when comparing the transitional 
quadrupole moment Qt [HI j even though the value of Q t 
for the TSD3 configuration is not much different. How- 
ever, if a TSD3 configuration is assigned, it appears to 
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correspond to an unrealistic increase of the spin values in 
the band compared with the values which appears most 
realistic from an experimental point of view [171 ] . The as- 
signment of a negative 7 configuration is in disagreement 
with Ref. [22I where it is concluded that a TSD2 mini- 
mum would be instable towards the TSD1 minimum in 
the tilted axis degree of freedom. This conclusion how- 
ever requires that the TSD1 minimum is (considerably) 
lower in energy than the TSD2 minimum which is the 
outcome of the CNS calculations as well as the calcu- 
lations of Ref. [22]. However, the relative energies of 
the three TSD minima are clearly uncertain within at 
least ±1 MeV. Thus, considering only the calculated en- 
ergies, it could be a configuration in any of the three TSD 
minima which should be assigned to observed yrast TSD 
band in 158 Er. One may also note that with present inter- 
pretations, there is no strong relation between the TSD 
bands in 158 Er and 168 Hf. 



VII. CONCLUSIONS 

In this study, we have made some modifications when 
solving the Hamiltonian in configuration-dependent 
cranked Nilsson-Strutinsky formalism in order to test the 
accuracy of some approximations. The nucleus 168 Hf is 
used as a test case and the observed highest-spin bands 
in this nucleus are analyzed. 

An important feature of the CNS formalism is that the 
off-shell matrix elements in the rotating harmonic oscil- 
lator basis are neglected in order to make it possible to 
fix configurations in more detail. For the yrast states, 
it is however straightforward to include all couplings and 
compare with the approximate CNS results. We conclude 
that the neglect of the off-shell elements of the Hamilto- 
nian matrix is acceptable; i.e. for the deformations and 
rotational frequencies which are reached in 168 Hf, it only 
leads to small errors which are essentially negligible com- 
pared with other uncertainties. This is also true for the 
cut-off error due to the limited number of oscillator shells 
in the basis, where our calculations show that if more 
than 9 shells (M m ax > 8) are included, both the total 
discrete energy and the total smoothed energy decrease, 
which means that the shell energy and thus also the total 
rotational energy is almost unchanged. 

The total energy was minimized in the full hexade- 
capole space indicating that for axially symmetric shape, 
it is generally sufficient to include only the standard £40 
degree of freedom which does not break the axial symme- 
try. On the other hand, at large triaxial deformation, it 
appears necessary to minimize the energy in all three hex- 
adecapole deformation parameters, (£40, £42, £44)- The 
study of the deformation space shows that in 168 Hf, the 
energy of the axially symmetric bands with normal defor- 
mation could as well be minimized in the restricted de- 
formation space (£ 2 , 7, £4), while the energy of the TSD 



bands should be minimized in the deformation space 
(£2, 7, £40, £42, £44)- We used these results and studied 
the structure of the experimental bands band 1, band 3, 
band ED and bands TSD1 and TSD2 in 168 Hf. 

In our studies of the high-spin bands of 168 Hf, the gen- 
eral conclusions are the same as in Refs. [H, EH but 
still with some important differences. Thus, we conclude 
that the crossing observed around / = 40 in the normal 
deformed bands are created when a particle is excited 
to the down-sloping 1/2[541] orbital but rather from the 
7/2[523] orbital and not from the 9/2[514] orbital sug- 
gested in Ref. 16]. The different conclusions are un- 
derstood from the fact that the number of particles is 
preserved in the unpaired formalism, while quasi-particle 
excitations with a fixed Fermi energy [l6| could lead to 
significant changes in the number of particles. More im- 
portant is however our suggestion that the ED band is 
built with two holes in the extruder hun orbitals. With 
such an excitation, the configuration is clearly different 
from that of the normal deformed bands with a calcu- 
lated transitional quadrupole moment in closer agree- 
ment with experiment. Because this configuration is not 
calculated as yrast, it is straightforward to analyze it 
only in formalisms where excited configurations with the 
same quantum numbers can be distinguished, e.g. by 
the number of particles (or holes) in high-j orbitals. For 
the TSD1 and TSD2 bands, we conclude that they are 
formed in a strongly deformed triaxial minimum with 
several particles excited to high-j intruder orbitals. This 
agrees with the assignment in Ref. [15| what concerns 
TSD1. However, while only the high-j particles were 
considered in that reference, the holes in the extruder or- 
bitals are as important according to our analyses, where 
it is mainly the neutron A^ = 4 orbitals which induces 
the triaxial shape. 

As a background for the study of the Hf bands, we 
did also investigate the filling of the orbitals of the TSD 
bands in Lu isotopes. These bands are naturally under- 
stood as having a proton configuration with one 243/2 or- 
bital occupied, compared with two such orbitals filled in 
the larger deformation TSD bands of 168 Hf. The neutron 
configurations in the N = 90 — 96 range are characterized 
by a successive filling of down-sloping orbitals in a region 
of low level-density which is created below four extruder 
orbitals emerging from the hxi/2 and JV = 4 shells, see 
Figs. Q] and 

The relative alignments, i e ff [33], between the TSD 
bands in the different nuclei was analyzed where the 
general features are understood from the contribution of 
the different orbitals which become occupied. Especially, 
it was concluded that the specific features of the yrast 
ultrahigh-spin band in 158 Er are best understood if this 
band is built in the TSD minimum with 7 < 0, i.e. for 
rotation around the intermediate axis. 
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